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Extracting the fundamental parameters 

Dirk Zerwas a 

a LAL, CNRS/IN2P3, B.P.34, 91898 Orsay Cedex, France 

If supersymmetry is discovered at the LHC, the extraction of the fundamental parameters will be a formidable 
task. In such a system where measurements depend on different combinations of the parameters in a highly 
correlated system, the identification of the true parameter set in an efficient way necessitates the development 
and use of sophisticated methods. A rigorous treatment of experimental and theoretical errors is necessary to 
determine the precision of the measurement of the fundamental parameters. The techniques developed for this 
endeavor can also be applied to similar problems such as the determination of the Higgs boson couplings at the 
LHC. 



1. INTRODUCTION 

In the following two scenarios for the extrac- 
tion of the fundamental parameters from experi- 
mental data will be discussed. The first example 
is vintage supersymmetry as an example for the 
precision achievable at the LHC when supersym- 
metry provides a multitude of signatures. The 
second example is difficult supersymmetry, be it 
that most, but not all of the supersymmetric part- 
ners are not observed, or that only a light Higgs 
boson is observed and its properties measured at 
the LHC. 

In a supersymmetric theory, each fermionic de- 
gree of freedom has a bosonic counter part and 
vice versa. Squarks, sleptons, neutralinos and 
charginos are the partners of the quarks, leptons, 
neutral and charged gauge and Higgs bosons. Su- 
persymmetric theories have "no" problems with 
radiative corrections, they predict a light Higgs 
boson with a mass less than 150 GeV, provide 
interesting phenomenology at the the TeV scale 
and may provide a link to the Planck scale. 

Many different models describe the phe- 
nomenology of supersymmetry. In the follow- 
ing sections, three of these will be discussed: 
mSUGRA as an example of a model with few pa- 
rameters, most of which are defined at the grand 
unification scale (GUT), Decoupled Scalars Su- 
persymmetry (DSS) as an example where a part 
of the supersymmetric spectrum is unobservable 



at the LHC, and the MSSM as a an example of a 
model defined at the clcctroweak scale with many 
parameters. Determining this model (in contrast 
to mSUGRA) may allow to measure grand uni- 
fication of the supersymmetric breaking parame- 
ters. 

R-parity is conserved in the following thus 
supersymmetric particles (cascade-)decay to the 
lightest supersymmetric particle (LSP). In the 
following the lightest neutralino will play the role 
of the LSP. It is stable, neutral and weakly inter- 
acting. The LSP is a candidate for dark matter. 
The experimental signature for supersymmetry at 
colliders is missing transverse energy due to the 
presence of the undetected LSP(s). 

In the present absence of supersymmetric sig- 
nals it is necessary to define benchmark points 
in order to study the potential of supersymmetry 
discovery and measurement. The reference point 
SPSla pQ has been studied intensely in the past 
years in LHC as well as ILC simulations [2]: In 
SPSla gluinos and squarks have masses of the or- 
der of 500-600 GeV, light sleptons have masses 
of 150 GeV and the Higgs boson is at the LEP 
limit. 



2. DETERMINATION OF PARAME- 
TERS 

The difficulties of determining the fundamental 
parameters from measurements can be illustrated 
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by taking mSUGRA as an example. The mass of 
the smuon depends on m , mjy 2 and tan/3. The 
mass of the chargino depends on va.x/2 an d tan /?, 
thus the same parameter has an impact on differ- 
ent measurements. Additionally the experimen- 
tal errors can be correlated and each theoretical 
prediction also has an error associated to it. Thus 
in order to disentangle the system to obtain the 
best possible precision on all parameters, a global 
Ansatz is necessary. 

For each observable in addition to a precise "ex- 
perimental" determination, a precise theory pre- 
diction must be associated. This means the most 
up to date calculations for masses, branching ra- 
tios, cross sections as well as dark matter predic- 
tions have to be used. A summary of different 
tools is listed in [3]. 

The first study of this type was performed 
in [4]. The Fittino [5] and SFitter [EE] collabo- 
rations studied the reconstruction of parameters 
with sophisticated techniques which will be de- 
scribed later. 

3. PREDICTIONS FROM PRESENT 
DATA 

While no direct observation of a supersymmet- 
ric particle is available, the theoretical predictions 
of precision observables are sensitive, via radia- 
tive corrections, to the parameters of supersym- 
metry. A wealth of precise measurements such 
as the mass of the W boson and the top quark 
mass have been made at LEP [8] and the TeVa- 
tron [9]. Additionally precise measurements in 
the b-sector, such as the branching ratio of b— >S7 
have been performed. The anomalous magnetic 
moment of the muon ((<? — 2) M ) is known with 
very good precision and WMAP has provided a 
measurement of the relic density. 

The comparison of the predictions with the 
measurements allows to determine the parame- 
ters of mSUGRA [ID]. Figure [T] shows the best- 
fit point in the (uq^/mo) plane. Due to (g — 2) M 
a positive sign of (i is favored. It is intriguing to 
note that the SPSla point, a point considered to 
be overly optimistic at the time of its definition, 
is close to the best-fit point. 

This approach can be taken a step further by 



without WMAP 




Figure f. The best fit result is shown in the 
(m^/mo) plane as well as the contours including 
and excluding the WMAP measurements. 



fitting the phenomenological MSSM, a weak scale 
model with about 20 supersymmetric parameters. 
Using a Bayesian approach with a linear prior and 
a logarithmic prior, one can infer on the MSSM 
parameters as shown in [llj . The dependence 
of the results on the prior should disappear as 
more information (such as supersymmetric mass 
measurements) become available. It is interest- 
ing to note that in this study the preference for a 
positive sign of the [i parameter from the muon 
anomalous magnetic moment is not observed. 

4. VINTAGE SUPERSYMMETRY 

The point SPSla provides a vast number of 
measurements at the LHC. In particular, the long 
cascade decay qx — » X2I ~ * £r&Q ~ * ^QXi can be 
observed well above the Standard Model and su- 
persymmetric backgrounds. The final state con- 
sists of opposite-sign same flavor leptons, i.e., 
electrons and muons, and hard jets. In this de- 
cay chain edges and thresholds can be measured 
by reconstructing invariant masses of different 
combinations: lepton-lepton, lepton-jet, lepton- 
lepton-jet. These endpoints are analytical func- 
tions of the masses of the particles, they do not 
depend on the underlying theoretical model. 

The results of the LHC study assuming an in- 
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Table 1 

LHC measurements [2] in SPSla are shown. The nominal values are calculated with SuSpect. Statistical 
errors, systematic errors from the lepton (LES), jet energy scale (JES) and theoretical errors arc all given 
in GeV. 
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tegrated luminosity of 300 fb _1 are shown in Ta- 
ble [TJ Several of the measurements are limited by 
the systematic error on the energy scale, which is 
assumed to be of the order of percent for jets and 
at the per mil level for leptons. Using the kine- 
matic formula for the edges and endpoints, the 
absolute masses of the supersymmetric particles 
can be determined. This procedure, usually per- 
formed with fits on toy experiments introduces 
additional correlations. For the determination of 
the fundamental parameters it is therefore prefer- 
able to start directly from the edges and thresh- 
olds. Simplifying the results of the analysis for 
the ILC, as soon as the particles are kinemati- 
cally accessible, they can be measured an order 
of magnitude more precisely than at the LHC. 

4.1. mSUGRA 

Two issues have to be addressed to determine 
the fundamental parameters from the experimen- 
tal measurements: finding the true/correct pa- 
rameter set from a strongly correlated system of 
measurements and determining accurately the er- 
rors on the parameters. 

For the first issue several different techniques 
have been developed to efficiently sample a multi 



dimensional parameter space. Fittino [5] has 
developed the technique of simulated annealing, 
which allows to cross potential boundaries. These 
boundaries could confine a simple fit based search 
to secondary minima. SFitter |6l7j has devel- 
oped weighted Markov chains which have the ad- 
vantage of efficient sampling in high dimensions. 
Markov chains are linear in the number of param- 
eters in contrast to a simple grid based approach 
where all parameters are sampled with a prede- 
fined step size. 

After having produced a full dimensional exclu- 
sive likelihood map, two types of projections are 
possible to reduce the number of parameters, e.g., 
to illustrate correlations between two parameters. 
The Bayesian approach introduces a measure but 
conserves the property of probability density in 
the projection. The frequentist approach uses a 
profile likelihood which ensures that the absolute 
minimum is always retained in the projection. 
In illustration of a two dimensional correlation 
plot is shown in Figure [2] for the frequentist ap- 
proach in mSUGRA. The correct parameter set 
was found from the measurements. 

The Markov chains also allow to identify sec- 
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Table 2 

Expected errors on the mSUGRA parameters at the LHC and ILC. Flat theory errors are used. 
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Figure 2. The result of the profile likelihood in 
the (m 1 / 2 /m.o) plane is shown [7]. The best-fit 
point is clearly identified and corresponds to the 
true SPSla parameter set. 



ondary minima. Even in a tightly constrained 
model such as mSUGRA secondary minima oc- 
cur. In particular, when in addition to the funda- 
mental parameters the Standard Model parame- 
ters are also allowed to vary within their exper- 
imental error. SFitter showed [7] that the inter- 
play of the top quark mass and the tri-linear cou- 
pling can lead to additional solutions. However 
these can be discarded easily by comparing the 
X 2 value of the solutions which are much larger 
than for the correct parameter set. 

To address the second issue, determination of 
the errors on the parameters, one can either de- 
termine them in a single fit using for example MI- 
NOS or by using toy-experiments, i.e., performing 
the parameter determination for a large number 
of datasets which have been smeared according 



to the experimental and theoretical errors. In the 
latter case the width (RMS or a Gaussian fit) of 
the distribution of the central value of the fits is 
the error. This procedure is more robust than a 
single fit and used most of the time, especially in 
the presence of correlations. 

While the experimental errors can usually be 
treated as Gaussian with well defined correla- 
tions, e.g., the energy scale error is fully corre- 
lated among measurements, the theoretical error 
deserves special attention. Here the RFit [12] 
approach is followed. Within the theoretical er- 
ror, the contribution to the % 2 is zero. Beyond 
this region the usual x 2 contribution calculated 
with the experimental error is used. The pro- 
cedure ensures that no particular value within 
the theoretical error region is privileged. The 
typical theoretical errors used for this study are 
3% for strongly interacting particles (squarks and 
gluinos) and 1% for weakly electromagnetically 
interacting particles (sleptons, neutralinos and 
charginos). In the Higgs sector an error of 2 GeV 
is used. 

The results are summarized in Table [2] for the 
LHC, ILC, the combination of LHC and ILC with 
and without theoretical errors. The precision of 
the LHC alone is at the level of percent for the 
determination of the parameters. It is improved 
by the ILC by about an order of magnitude. In- 
cluding the theoretical errors has an impact on 
the precision at both machines, the errors are 
larger by a factor of three to four. Thus the preci- 
sion of the parameter determination at the LHC 
is limited by the precision of the theoretical pre- 
dictions. The SPA project [13] aims to improve 
the situation. It is also important to note that at 
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the LHC the precision on the top quark mass of 
1 GeV has an impact of about 10% on the preci- 
sion of the determination of the parameters. 

4.2. MSSM 



0.006 




400 600 800 1000 
M 1 

Figure 3. The result of the profile likelihood 
for the gaugino mass parameter Mi is shown [7]. 
Several peaks are identified corresponding to sec- 
ondary maxima. 



The search for the global minimum is more 
complex in the MSSM where more parameters 
have to be determined. If the unification of the 
first and generation sfermions is not assumed, 
19 parameters have to be determined. A com- 
plex mix of different techniques has to be applied: 
markov chains and a gradient search to refine the 
resolution of the minima. The result of the pro- 
file likelihood procedure at the LHC is shown in 
Figure for the gaugino mass parameter Mi . 

Several secondary minima are observed. They 
cannot be distinguished by their \ 2 value which 
is identical in all cases. The degeneracy of the 
eight secondary minima can be understood qual- 
itatively by analyzing the gaugino-Higgsino sec- 
tor. At the LHC in SPSla only three neutrali- 
nos can be measured according to the studies 
available today, but the sector is determined by 
four parameters Mi, M 2 , (i and tan/3, neglect- 
ing radiative corrections. The system is there- 
fore under-constrained. Additionally Mi and M 2 
can be exchanged without a change in the LHC 



observables. The situation is improved dramat- 
ically when the ILC measurements are added as 
they complete the neutralino sector and add the 
masses of the charginos providing additional con- 
straints to determine the MSSM parameters with- 
out ambiguity. 

The MSSM is defined at the electroweak scale. 
The definition of its parameters does not depend 
on the model of supersymmetry breaking. Thus, 
having determined the parameters at the weak 
scale, the extrapolation of the parameters to the 
GUT scale can be performed. Instead of assuming 
unification at the GUT scale as in mSUGRA, it 
can now be tested. As eight degenerate solutions 
have been determined the extrapolation is per- 
formed for each one of them. Two examples are 
shown in Figure |4] (top) . Figure H (top left) cor- 
responds to a solution which is the correct SPSla 
MSSM parameter set. GUT unification of the 
gaugino masses is observed as expected. How- 
ever in Figure [4] (top right) GUT unification is 
not observed. Of the eight degenerate solutions 
only one unifies the gaugino masses at the GUT 
scale, one almost does (this corresponds to the 
parameter set where only the sign of \x is changed 
with respect to SPSla) and the six other solutions 
clearly do not unify. 

In addition to the gaugino sector the extrapo- 
lation of the scalar sector in the first two genera- 
tions can also be studied at the LHC. The third 
generation is under-constrained due to the lack 
of measurements in the stop sector. In Figure [4] 
(bottom) the extrapolation is shown on the left 
for the correct SPSla MSSM parameter set and 
on the right for the same set with the exception 
of the stau parameter which has been moved far 
from the true solution (the \ 2 values are identi- 
cal). Unification can be observed in the first case, 
but not in the second one, while the gaugino uni- 
fication is observed in both cases. This can be un- 
derstood from the RGE equations. The RGEs in 
the gaugino sector essentially decouple, whereas 
in the RGEs of the scalar sector for a parameter 
all other parameters. A bad measurement, i.e., 
far from its true value, therefore has a large im- 
pact on the extrapolation of all scalar parameters. 
The extrapolation is stabilized by the measure- 
ments of the ILC. The combination of LHC and 
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Figure 4. Top: Extrapolation of the inverse of the gaugino mass parameters to the GUT scale for two 
of the degenerate solutions at the LHC. Grand unification is observed in one case corresponding to the 
correct parameter set (left), but not in the other solution (right). Mi is in black, M2 red and M3 is in 
green. Bottom: Extrapolation of the first and second generation scalar mass parameters to the GUT 
scale for two of the degenerate solutions. Grand unification is observed in one case (left) corresponding 
to the correct SPSla MSSM parameter set, but not the other one (right). 



ILC therefore will allow to measure grand unifi- 
cation of all supersymmetric parameters, whereas 
for the LHC, without further measurements, one 
can only observe that one of the ambiguous solu- 
tions is compatible with GUT unification. 

4.3. RELIC DENSITY 

At this point all supersymmetric parameters 
have been determined, therefore the complete 
particle spectrum can be deduced. From the spec- 
trum and its couplings the relic density is pre- 
dicted [M]. Neglecting the theoretical error and 
performing the analysis in mSUGRA, a precision 
of about 2% is expected at the LHC with an im- 



provement by an order of magnitude expected 
from the ILC. The precision depends on the exact 
nature of the parameter set and is not valid in all 
scenarios. The precision compares well to the one 
expected by the Planck satellite which was suc- 
cessfully launched recently. Thus the confronta- 
tion of the collider prediction of the relic density 
and the measurement of the cosmic microwave 
background fluctuation will provide for interest- 
ing physics studies in the future. 
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Table 3 

Signatures included in the Higgs coupling analysis for a Higgs mass of 120 GeV and an integrated 
luminosity of 30 fb _1 . The factor after the background rates describes how many events are used to 
extrapolate into the signal region. The last two columns give the one-sigma experimental and theory 
error bars on the signal. 
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5. DIFFICULT SUPERSYMMETRY 

While one can hope for the LHC discovery of 
SPSla-like supersymmetry as hinted by the anal- 
ysis of precision observables of the electroweak 
sector, b-sector and the relic density, difficult sce- 
narios must also be envisaged. One of these is 
DSS, also known as Split-SUSY ( [16117] and refer- 
ences therein) . In these models the scalars are too 
heavy to be produced at the LHC. However sev- 
eral observables are still present: the mass of the 
lightest Higgs boson, the mass difference between 
the second lightest neutralino and the LSP, the 
cross section of the tri-lepton signal, ratios of the 
decay of supersymmetric particles to the Z boson 
and the gluino production cross section. Putting 
all of this together, the DSS parameters can be 
determined and measured at the LHC. This is an 
indication that the LHC can also handle difficult 
scenarios. 

5.1. HIGGS 

Another difficult scenario at the LHC would 
be that only a light boson with a mass of about 
120 GeV is discovered at the LHC but no (other) 
new particles. The expected measurement chan- 



nels and their experimental and theoretical errors 
are shown in Table [3] While an ILC is capable 
of a model independent search for the Higgs bo- 
son, at the LHC the search and measurement is 
restricted to well defined final states. Studying 
the precision of the determination of the Higgs 
couplings could shed light on new physics. 

The study presented here [19] uses essentially 
the same MonteCarloData as 18 with two im- 
portant differences: In agrement with more recent 
experimental studies by ATLAS and CMS [20121] . 
the ttH— »bb signal is reduced by 50%. On the 
other hand the new theoretical analysis of the 
ZH(— >bb) (subjet analysis) [22 is used, confirmed 
by ATLAS within 10% [23]. 

The theoretical errors on the production cross 
section are between 7% and 13%. The errors on 
the branching ratio calculations are between 1% 
and 4% (c quarks). The experimental errors have 
a statistical uncorrelated component and a sev- 
eral sources of systematic error which are highly 
correlated among measurements. Determining 
the Higgs couplings from this system is therefore 
similar to the determination of supersymmetric 
parameters, i.e., the same techniques can be used 
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Figure 5. The result of the profile likelihood is shown for the Htt versus the HWW couplings on the (left) 
without allowing additional contributions to effective couplings and on the (right) allowing additional 
contributions to effective couplings. 



to extract the couplings. The additional difficulty 
to be mastered as some of the channels are in 
the Poisson regime is the convolution of the flat, 
Gaussian and Poisson errors. 

The parameters are defined in the following 
way for tree level couplings: 

9hh — </;,// (i + A JjH ) (i) 

as deviations from the Standard Model value. For 
the loop induced couplings such as the 77H and 
ggH couplings, the definition is as follows 

Hun -/^ (1 • • A„„) (2) 

The additional term A^j^ is the modification of 
the jjH effective coupling induced by a deviation 
of the tree level couplings of the particles in the 
loop. 

As in the supersymmetric case, first a full 
dimensional likelihood map is calculated from 
which the Bayesian and frequentist projections 
are possible. In the Higgs analysis however no 
true secondary minima exist and the noise effects 
of the Bayesian integration are large, essentially 
washing out the signal. Therefore only the fre- 
quentist projection is pursued. 

Two main observations can be made after the 
projections to two dimensions: there is a general 
positive correlation among all couplings with the 
exception of the bbH coupling. This is due to the 
fact that for a Higgs boson mass of 120 GeV the 



total width is dominated by the decay to b-quarks 
(approximately 90%). The total width enters in 
the denominator of every observable, whereas the 
individual partial widths, proportional to the cou- 
plings (squared) enter in the numerator. 

The second observation, as shown in Figure [5] 
(left), is that as long as no additional contribu- 
tions in the loop induced couplings are allowed 
(AjjH Axed to 0), a preference for the correct 
sign of the ttH coupling with respect to the WWH 
coupling is observed. This is due to the Higgs de- 
cay channel to photons. Allowing genuine anoma- 
lous couplings in the loop induced couplings de- 
stroys this distinction completely as the addi- 
tional parameter allows for the compensation of 
the sign preference without a penality to be paid 
in the % 2 value. 

It is also illustrative to remember what the 
LHC cannot do in addition to what it can do. All 
observables are have the structure (simplified) of 
cross section times branching ratio. Expressed in 
partial widths we have TnH^jjH /^tot- As partial 
widths are proportional to the square of the cou- 
plings, the structure becomes g 2 g 2 /T tot . Thus if 
a coupling not measured at the LHC, e.g. Hcc is 
much larger than expected, the total width will 
increase, but the coupling determination would 
show deviations from the Standard Model, i.e., 
couplings smaller than expected. A simultane- 
ous determination of the individual couplings and 
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Figure 6. The impact of the subjet analysis on the Higgs boson coupling determination is shown. On the 
(left) the full sensitivity is used, on the (right) the channel is removed completely. 



the total width is not possible as no model in- 
dependent channel is available. Thus the LHC 
can measure absolute couplings only under the 
assumption that all non-measured couplings are 
equal to their Standard Model prediction. 

With this caveat in mind, assuming an inte- 
grated luminosity of 30 fb _1 , the Higgs boson cou- 
plings can be determined with a precision of 30% 
(WWH) to 50% (ttH). As before toy experiments 
were used to obtain these results. The precision is 
somewhat affected [19] by whether effective cou- 
plings are allowed to vary or not. The effect is 
of the order of up to 20% (relative) . If the ratios 
of couplings are analysed instead of the couplings 
and the WWH coupling is used as reference, the 
precision is of the order of 30%. While part of 
the theoretical and experimental error cancel and 
therefore one could expect some improvement, 
with an integrated luminosity of 30 fb _1 the sta- 
tistical error is still dominant. The conclusions 
could therefore change as more data will become 
available. 

In Figure M (left) the (bbH,WWH) coupling 
plane is shown using the full sensitivity of the 
analysis. In the same Figure (right) the impact of 
removing the analysis completely is shown. The 
spread of the central region is much larger. The 
subjet analysis is essential for the Higgs coupling 
determination at the LHC because of the impor- 
tance of the bbH coupling for low Higgs boson 



masses. 

Once the Higgs couplings are determined one 
can ask whether the precision is sufficient to ex- 
clude new physics. The task looks daunting given 
the precision of the individual couplings. How- 
ever not only the individual couplings are impor- 
tant but also their correlations. The loglikelihood 
is used as an estimator take all correlations into 
account. 

In the gluophobic Higgs [24] model the stop 
quark and top quark contributions cancel to re- 
duce the ggH coupling to 24% of the Standard 
Model value, i.e. the cross section for gluon fu- 
sion processes is reduced by a factor 25. In this 
case at 90% C.L. 46% of the toy experiments are 
not described by the Standard Model. 

In a second scenario the parameter set SPSla 
was moved out of the decoupling region by mod- 
ifying the mass of the CP-odd Higgs boson A, 
the top tri-linear coupling and tan/3. Some of 
the channels are enhanced, others are reduced by 
almost 40%. At 90% C.L. 77% of the toy experi- 
ments are not described by the Standard Model. 

6. CONCLUSIONS 

If vintage supersymmetry, as hinted by the 
electroweak fits, is found at the LHC, a large num- 
ber of measurements will be available which can 
be translated into the fundamental parameters of 
mSUGRA and even the MSSM. To unambigu- 
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ously determine the parameter set and measure 
grand unification of the breaking parameters how- 
ever the ILC will be necessary. The LHC is pre- 
pared to determine the fundamental parameters 
in difficult scenarios, e.g., if only a Higgs boson 
is discovered at the LHC. A determination of the 
Higgs couplings can shed light on new physics. 

The extraction of the fundamental parameters 
of any theoretical model is a formidable task 
which requires a close collaboration of experimen- 
talists and theorists to develop the sophisticated 
analyses to determine the fundamental parame- 
ters at the LHC (and beyond at the ILC) what- 
ever nature has in store. 
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